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Fluctuations, correlations and non-extensivity 
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The present status of investigations on fluctuations and correlations seen in high energy multiparticle produc- 
tion processes made using the notion of nonextensivity is reviewed. 
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I. INTRODUCTION 

There are many examples that distributions of produced 
particles deviate from the expected exponential form towards 
its power-like generalization [VJ, called the ^-exponential and 
know also as Tsallis distribution ]2[], 



be described by gamma distribution, 
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described by parameter q (in cases of interest to us X = 
pr or X = mrcoshy). Distributions of this type were pro- 
posed on phenomenological grounds long time ago as useful 
parametrization of data interpolating between results coming 
from soft and hard scattering [3]. On the other hand, when 
seen from the information-theory point of view or from the 
thermodynamical approach, such distributions arise from the 
use of non-extensive Tsallis entropy [4] (indexed by q) instead 
of the usual Shannon-Gibbs-Boltzmann one: 
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SsBG = -YtPi^Pi- ( 2 ) 



Notice that S q is nonextensive in such sense that for two in- 
dependent systems, A and B (in the usual sense, i.e., that 
the corresponding probabilities of their occurrence factorize, 
pij(A+B) = pi{A)pj(B)) one has that 

S q (A+B) = S q (A)+S q (B)+(l-q)S q (A)S q {B). (3) 

The entropic index q is thus also a measure of nonextensivity 
in system under consideration and therefore is also called non- 
extensivity parameter. 

What are the circumstances leading to q ^ 1 ? The common 
consensus is that it happens whenever [4]: - there are long 
range correlations in the system (or system is small - notice 
that our Universe is small with respect to the gravitational in- 
teractions!); - there are memory effects of different kind; - the 
phase-space in which system operates is limited or has fractal 
structure and, finally, - there are intrinsic fluctuations in the 
system under consideration. 



II. THE ROLE OF INTRINSIC FLUCTUATIONS 

Our works, which I shall review now, were concerned with 
the last possibility mentioned above, namely that 1 1 — q\ mea- 
sures fluctuations in the hadronizing system under considera- 
tion. As shown in [5] in the case when these fluctuations can 
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where (...) are the respective averages taken with respect to 
/. In general case one refers to the concept of the so called 
superstatistics discussed at length in_|0]. This approach has 
been recently further generalized in |7|] . 

The situation encountered can be visualized in the follow- 
ing way. Already in |8J] it was suggested that the overwhelm- 
ingly success of thermodynamical models has its origin in the 
fact that out of the large number N of secondaries produced 
in a given event only one (very rarely two) is chosen for mak- 
ing distributions. Inevitably, the remain (N — 1) act as a kind 
of heat bath which action can in most cases be described by 
single parameter T (called usually temperature because of as- 
sociations with thermal models). In this case one has simple 
exponential distributions (i.e., q = 1). However, because such 
"heat bath" is in reality neither homogeneous nor infinite, it is 
natural to expect that it should be described by more param- 
eters. The simplest extension is to assume that T fluctuates, 
replace T — > To = (T) and introduce another (second) param- 
eter describing its fluctuations - this leads to q ^ 1 in most 
natural way. Mathematically it can be realized by introduc- 
ing to the well-known linear Langevin equation with additive 
noise, modelling the Brownian motion of some test particle 
(and leading to a Boltzmann distribution as stationary solution 
of the corresponding Fokker-Planck (FP) equation), a small 
multiplicative noise. This leads immediately to Tsallis dis- 
tributions of relevant observable as exact stationary solution 
of the corresponding FP equation with (q— 1) being directly 
connected with the strength of this new term (cf. ||5l |7[]). To 
close this point it must be mentioned that the role of corre- 
lations is so far not excluded but it seems that in the realm 
of high energy multiparticle production they are, at least at 
circumstances considered so far, not the dominant component 

Experimental data of interest come usually in two cate- 
gories: as distributions in longitudinal phase space (usually in 
rapidity y) averaged over transverse momenta, dN/dy, or dis- 
tributions in transverse momenta, dN/dpj, averaged over y or 
taken for some window in y. We have therefore also two kinds 
of "temperatures", Ti and Tj, and their fluctuations are thus 
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described by, respectively, q^ and qj. Out of these two classes 
only exponential shape of dN/dpr spectra is usually regarded 
as indication of thermal-like character of such processes with 
a kind of local thermal equilibrium setting in and character- 
ized by temperature T — Tj, which can be deduced from the 
corresponding slopes of dN/dpr- Any deviation from them 
is then interpreted as signal for some dynamical (nonequilib- 
rium) effects showing up (like, for example, the flow or decay 
of resonances, see iflltl and references therein). Instead of try- 
ing to exclude them one can investigate the possibility that the 
observed nonexponential spectra result from some new form 
of equilibrium characteristic of nonextensive thermodynam- 
ics (incorporating, among others, also effects of the aforemen- 
tioned factors). The plausibility of such approach is illustrated 
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FIG. 1: Examples of Tsallis distributions: (a) fit to rapidity spectra 
for charged pions produced in pp and pp collisions at different ener- 
gies, (b) fit to pj spectra from UA1 experiment (see 11211 for details 
and references to experimental data). 

in Fig. Q] where examples of both types of distributions are 
shown 1 12]. The characteristic feature encountered here is that 
1l > It and that Tj, > > Tj {qt ranges in linear fashion from 
1.05 for y/s = 20 GeV to 1.33 for y/s = 1800 GeV whereas 
variation of qj is limited to between 1 .095 for y/s = 200 GeV 
and 1.11 for a/s = 900 GeV; the respective changes of Tl are 
from 1.76 GeV to 55.69 GeV whereas those of Tj are much 
smaller varying from 0.134 GeV to 0.14 GeV). The immedi- 
ate question arises: what is the meaning of qi, or, equivalently, 
what is the meaning of fluctuations of "partition temperature" 
T\P. The answer is that qi = l + l/k, where k is parameter 
defining (in addition to mean multiplicity («)) the Negative 
Binomial distributions of multiplicity observed in such col- 
lisions 1 16]. It is because Tl ~ M/(n), where M is energy 
available for production of secondaries, therefore for M kept 
constant fluctuations of Tl mean automatically fluctuations of 
(n) and those lead immediately to NBD 01611 . Interestingly 
enough NBD arises also in a natural way when particles are 
produced with energies distributed according to Tsallis dis- 
tribution (when they are distributed according to Boltzmann 
distributions one gets Poisson distribution instead, see 191]). 

Let us proceed now to recent RHIC data on pj distribu- 
tions, which consist our potential source of the parameter T, 
as mentioned above. As shown in 11311 . using approach based 
either on nonextensive statistics or on stochastic ideas one 
can successfully account for the whole range of the observed 



transverse momenta. This is because in both cases the resul- 
tant distributions are intrinsically non-exponential. In 01411 the 
same data were analyzed using Hagedorn model with spec- 
trum of resonances given by p(m): 
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As seen in Fig. [21 although p(m) introduces already some 
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FIG. 2: Example of analysis of STAR data IT^l by using usual Hage- 
dorn formula (eq. ((7} with q = 1, left panel) and its nonextensive 
generalization (eq. with q > 1, right panel), see fl4ll for details. 

fluctuations to the system, experimental data can be fitted only 
with qj > 1 (of the order of qj ~ 1 .00015). The new and po- 
tentially very interesting fact is that similarly good fit can be 
obtained with simple ^-exponential (i.e., by putting p(ni) = 1 
in eq. (0), however, in this case q is noticeably greater and 
equal to qj = 1.065. This effect is seen for all RHIC data 
and for all centralities. It strongly suggests that including 
some well identified sources of fluctuations to simple statis- 
tical model (represented here by function p(m) as given by 
eq. (JHJ) accounts for some fluctuations and lowers therefore 
the value of parameter qj. However, as seen in Fig. |2] data are 
very sensitive to q and, as has been shown in IU4I1 . to obtain 
good estimation of To = 1 /po one has to resort to nonextensive 
version of eq. (O with qr > 1. 




FIG. 3: Example of EGE fnil: I nelasticity distributions %(K) (nor- 
malized to unity) obtained in 1 16] from analysis of multiparticle pro- 
duction data for *Js = 200 GeV and ^/s = 900 GeV fitted by gaussian 
(full lines) and lorentzian (dashed lines) distributions, respectively 
(see 1 16] for more details). 

We would like to close this Section by mentioning that there 
is yet another source of fluctuation in the production process, 
which also leads to nonextensivity but not of Tsallis type, i.e., 
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with gaussian distributions instead of eq. ((T). Their source 
lies in the fact that usually only fraction K (called inelastic- 
ity) of the total energy of reaction is used for production of 
secondaries and this fraction can fluctuate. It turns out that 
such fluctuations are most naturally described in the frame of 
the so called Extended Gaussian Ensemble (EGE) by simple 
gaussian distribution, %(K) ~ exp [—(K— (K)) 2 /(2a 2 )] , (cf. 
IU7I1 for details and references) and that they can reasonably 
fit distributions of %(K) for different energies obtained in Illql . 
see Fig. [3] In fact, closer inspection of Fig. [3]shows that still 
better fit can be obtained when using the so called lorentzian 
distribution, which is nothing else as g-gaussian with q = 2. 
It would mean then that even here there are still some other 
fluctuations seen in data, which are undisclosed by the usual 
EGE approach. 

III. SUMMARY 

To summarize: we argue (since quite a time already) that 
whenever one finds in data a hint for eq. (fTJ the immediate 
suspicion should be that they hide some intrinsic fluctuations 
with strength given by q — 1 . What is the dynamical origin of 
these fluctuations is, however, out of the scope of the proce- 
dure discussed here as q accounts for all of them. The q > 1 
results should then urge experimentalists to devise some spe- 
cial measurements devoted to search for fluctuations. 

To allow the reader to make personal judgement in what 
concerns the possible role of correlations let us list recent at- 



tempts to interpret parameter q by some dynamical correla- 
tions caused either by incomplete phase space occupance HH 
or by some specific changes introduced in generalized form of 
the Boltzmann equation (either by using corresponding colli- 
sion rates nonlinear in the one-particle densities or by using 
nontrivial energy composition rules in the energy conserva- 
tion constraint part, cf. |191 and references therein). It should 
also be added at this point that, as advocated in Il20ll . one can 
also view q > 1 as a general (leading order) finite-size effect. 

Finally, one should be aware of the fact that there is still 
an ongoing discussion on the meaning of the temperature in 
nonextensive systems. However, the small values of param- 
eter qj deduced from data allow us to argue that, to first ap- 
proximation, 7b can be regarded as the hadronizing tempera- 
ture in such system. One must only remember that in general 
what we study here is not so much the state of equilibrium but 
rather some kind of stationary state (see 12111 and references 
therein). 

Let us close by saying that this subject is still an open issue 
for further research (like, for example, event-by-event analysis 
of data l22ll or hydrodynamical models ll23ll ). 
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